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We consider the low-temperature limit of the long-distance asymptotic behav- 
ior of the finite temperature density-density correlation function in the one- 
dimensional Bose gas derived recently in the algebraic Bethe Ansatz framework. 
' Our results confirm the predictions based on the Luttinger liquid and confor- 

mal field theory approaches. We also demonstrate that the amplitudes arising 
in this asymptotic expansion at low-temperature coincide with the amplitudes 
associated with the so-called critical form factors. 
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1 Introduction 



The model of one-dimensional bosons interacting through a two-body <5-function potential is 
described by the Quantum Nonlinear Schrodinger equation (QNLS model). Recently [TJ, the 
Lagrange series method was applied in the framework of the algebraic Bethe ansatz to the 
derivation of the long-distance asymptotic expansion of the finite temperature density-density 
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correlation function in this model. In the present paper, we consider the low-temperature limit 
of this result. 

The QNLS model is given by the Hamiltonian 

L 

H= [ (d x ^d x ^ + c^^^$> - h^y) dx. (1.1) 



Here \E f and are Bose fields possessing canonical equal-time commutation relations, c > is 
a coupling constant and h > the chemical potential. The results obtained in [1] are relative to 
the thermodynamic limit L — > oo of this model. 

The density operator j(x) = ty'(x) ^(x) defines the operator of the number of particles in 
the interval [0,x] 

X 

Q x = j j(z)dz. (1.2) 
o 

The generating function for the density-density correlations reads 

„ T lim <5Mf^M, (1 . 3) 

were T is the temperature, a a complex number and \Qt) any eigenstate of H that goes to 
the state of thermal equilibrium in the infinite volume limit. Indeed, the correlation function of 
densities (j(x)j(Q))T can be obtained from (|1.3p as 

8tt 2 dx 2 da 2 



(j(-)j(0))T=-^^^(e 2 ^)T 



11.4) 



a=0 

We have shown in [1] that the large- a; asymptotic expansion of the generating function (jl .3|) 
(and respectively of the two-point function (jl.4p ) is given in terms of solutions to a set of non- 
linear integral equations closely related to ones arising in the quantum transfer matrix approach 
[2 El HI El El E] • Below, we solve these equations in the low-temperature limit. This computation 
allows us to reach two goals. 

On the one hand, one can argue that the Luttinger liquid [8] and conformal field theory 
(CFT) approaches [9| \10\ \TT\ [12j [131 El ESI EE] can be used to predict the large- a; asymptotic 
behavior of the low-temperature correlation functions in massless one dimensional quantum 
models. The QNLS model does belong to this class. Thus, our results give us a possibility to 
confirm these predictions by a direct derivation based on the algebraic Bethe ansatz. Namely, we 
show that in the low-temperature limit (T — > 0) the asymptotic expansion (x —> oo, xT —> oo) 
of the density-density correlation function takes the following form (at leading order for each 
oscillatory term): 

<j(aO;(0)>T = D 2 — +y Ai e 2 ^ ( gp-, . (1.5) 

UV ' K 11 2 S inh 2 (^/ Wo ) ,t=i \smh(7TTx/v )J V ' 
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Above appear several constants that will be defined in section [3l namely, the average density of 
the gas D, the Fermi momentum k F , the velocity of sound on the Fermi boundary vq and the 
value of the dressed charge at Fermi boundary Z. The coefficients Ai do not depend on T. All 
the dependence of the amplitudes on T has been gathered in the pre-factor {-kT / vq) 2 ^ 2 2,2 . 

This form is in full agreement with the CFT predictions. Moreover it provides one with an 
asymptotic behavior that is also valid in the full scaling region of xT and in particular in the 
T = case. 

On the other hand our approach allows us to calculate the constant coefficients Ag in ([1.5p . 
We show that Ai are related to the amplitudes of the so-called critical form factors intro- 
duced in |17| and arising in the study of the model at T = 0. More precisely, the coefficients 
A^{-kT / vq) 2(2 2,2 determined for the system in the thermodynamic limit and at small but finite 
temperature T are equal to the amplitudes of the critical form factors corresponding to umklapp- 
type excited states of momentum 2k F t and determined for the system of large but finite size L 
at zero temperature, with the identification vq/T i— > iL. We will show this coincidence by means 
of straightforward calculations carried in the core of this paper. 

This article is organized as follows. In section [2] we recall the results obtained in [I]. In 
section [3] we present the thermodynamics of the QNLS model at low temperature. In particular 
we solve the non-linear integral equation determining the asymptotic expansion of the correlation 
function in the low-temperature approximation. This allows us to obtain the rates of exponential 
decays in section U] and the constant amplitudes in section [5] The expansion (|1.5p is derived in 
section [6l Various estimates of the low-temperature behavior of the integrals that we deal with 
are gathered in three appendices. 



2 Long distance asymptotic behavior at general temperature 

The state of the thermal equilibrium in the QNLS model is described by the Yang- Yang equation 
|18| for the excitation energy e(A) 

e ( A ) = a 2 - h - £ J K(\ - n) log (l + e"^) dp, (2.1) 

R 

and the integral equation for the total density pt(\) 

Pt(X) - ^ J K{\ - iA)<&{ii)pt(p) dp = i-. (2.2) 

R 

The kernel K(\) and the Fermi weight $(A) appearing above read 

K ( A ) = A^' m = (l + e^y 1 - (2-3) 
Below the poles of the Fermi weight will play an important role. Therefore we introduce the 
roots {r*} of the equation 1 + e _£ ^ ^ T = 0, where rt (resp. rj) belong to the upper (resp. 
lower) half-plane (see Fig. [TJ. 
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The asymptotic expansion of the generating function (e 27ria ® x )y is given in terms of the 
solutions to non- linear integral equations similar to (|2.ip . Let us choose n points (n = 0, 1, . . . ) 
sj~ in the upper half-plane and n points sj in the lower half-plane. We then introduce non-linear 
integral equation for a function u(A) 

u(X) = A 2 - h a - ^ fK(X-ij)]og(l + e-^dii + iT^(9(X-^)-e(X-8j)), (2.4) 



where h a = h + 2iriaT and 



0(A) = t log (^±^1, 9'(X) = K(X). (2.5) 



ic — X J 

Below we will show that, in the low-temperature limit and for {s*} fixed, the solution to the 
equation (|2.4p always exists. Clearly, this solution depends on the parameters {s*}: u(X) = 
u(X\{s + }, {s~}). By imposing the constraints 

/ n(s±|{s + },{r})\ 

l + expf- jlt J' 1 J M =0, j = l,...,n, (2.6) 

we obtain a system of equations, which fixes the sets {s*}^ that are relevant for the description 
of the long-distance asymptotics. The subscript i enumerates these sets. The long-distance 
asymptotic expansion for the generating function (^e 27Tia ^ x )x can then be organized into a sum 
parameterized by the functions tij(A) = Ui(X\{s + }i,{s~}i) that solve (|2.4p with a corresponding 
set of roots {s ± }i: 

( e 2«faa.) T = J2e- xp[ui] B[ui] + o ( e - imaxS (pK])J , (2.7) 

i 

where and are functionals of Ui(X) whose explicit form will be specified later. 

Observe that equation (j2.4|) can be recast in the form 

u i Cm) \ 

dfi, (2.8) 



Ui (X) = X 2 -h a -^- j K(X-fx) log (l + e 



where the contour Cj is such that the roots {s ± } are located between the real axis and Cj. We 
also demand that the contours Ci separate the sets {s ± }i from all over possible roots of the 
equation 1 + e~ Wi ^' T = and from all the roots {r - *} (see Fig. [T]). Then one can interpret the 
asymptotic expansion of l^ 2,mOL Q- x ^ T as being given by the sum over the different possible choices 
of contours C{. 

The explicit expressions for the rates of exponential decay p[u{] and the amplitudes B[u{] 
can be written down in terms of integrals over the contours Cj. Let us introduce an auxiliary 
functions Zi(X) by 
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Figure 1: The roots {r*} are depicted by (o), the roots {s ± }i are depicted by (•). Other roots 
of the equation 1 + e~ Ui ^ s± ^ T = are depicted by (x). The contour Ci bypasses the roots {s + }i 
from above and the roots from below. It also separates the points {s^i from the other 

points {s ± } as well as from all the poles of the Fermi weight {r*}. 



Then the rates of exponential decays take the very simple form 

p[u k ] = i [ z k (X)dX. (2.10) 



The expressions for the amplitudes B[ui] are much more cumbersome. We present them in 
the form 

B[ui] = B d [ui]B.[ui], where B d [ Ui ] = exp ( J d\ df)j , (2.11) 

Ci 

and the symbol /x + means that the variable [i is slightly shifted to the left of the oriented 
integration contour Cj. The reason we separate the coefficients B\uj\ into two factors Bd[ui] and 
B s [ui] is that these go, in the low-temperature limit, to the discrete and to the smooth parts of 
critical form factors respectively, as defined in |17| . 

In order to describe -B s [itj] we first introduce the Cauchy transform operator on the 
contour Ci 

L 4 [*](«) = J J^dX, (2.12) 

and a functional Cq = Co[^j] 



Zi(X)zi(n) 
(A — fj, — icY z 



co&]= I ,r rv A2 d\d». (2.i3) 



Ci 
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Then 

det (i + it/W [zi]) det(l + [ Zi ] 



5 s M = (e 2 ™*-l)^ 



det [l det [/-!]((«.)] 



Lg , [^j] (#i +zc) 27via-\-L^ [zj\ (0\ —ic) 



-i 



^— Lg [zi](02— ic) [zj](02+ic) 



-1 



(2.14) 



The first line of (|2.14p contains a ratio of Fredholm determinants. The integral operators / — 
and / — -^K^ have the kernels 

1 + eT 1 + e t 

The first of these integral operators acts on the real axis and the second one acts on the contour 
The operators 1 + 2^ (1) M and J + 5H^ (2) N both act 

on a anticlockwise oriented closed 

contour surrounding Cj. Their kernels are given by 



E^W, M) = -e^ W(w) • .^""^"^jfV; ., , (2.16) 



and 



#< 2 W, N) = e -^N(»0 . ^^-.^-^(7-^) , v (2.17) 

v ' ' L lj; — L,s Zj (ui'-ic) 2-Kia-Lfl \z^\(w'+^c) , v y 

*«(A) = 3— - — - . (2.18) 
A + ic \ — ic 

Finally, observe that the kernels U^ 1,2 \w, w', [zi]) as well as the coefficient in the second line of 
(|2.14p depend on arbitrary complex numbers 0\ and O2 located inside of the contour where the 
operators U^' 2 ^(w, w' , [zi]) act. One can prove (see |19| ) that the total combination (|2.14p does 
not depend on the specific choice of these parameters. 



where 



3 Thermodynamics at low temperature 

In this section we provide a list of necessary formulae describing the thermodynamics of the 
QNLS model at low temperature. We begin our discussion with the T = case. 



3.1 Zero temperature 

It is known |18] that for a positive chemical potential h > the solution e(X) to the equation 
(|2.ip has two roots ±q(T) on the real axis: e(±q(T)) = 0. Hereby e(X) > for |A| > q(T) and 
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e(A) < for |A| < q(T). Let e(A) e (A) and q(T) -»• g as T -»• 0. Then, 



/ £( a)n f 0, |A|>g, 

limTlog 1 + e-— ) = { . (3.1) 

V ' { -eo(A), |A|<g, 

It is then straightforward to check that equation (|2.ip turns into a linear integral equation for 
the dressed energy £o(X): 

9 

£o(A) - ^ K(X - fj,)s Q (n) dn = X 2 - h, e (±q) = 0. (3.2) 



-9 



At T = the state of thermal equilibrium goes to the ground state of the QNLS model. The 
Fermi weight i?(A) (|2,3p turns into the characteristic function of the interval [—q,q]. Therefore 
the equation (|2,2p for the density takes the form 

9 

PtW -^Jk(\- M )ptOi) d// = ^, T = 0. (3.3) 

-9 

Another important characteristic of the ground state is the dressed charge Z(\). In the 
QNLS model it is proportional to the density Z{X) = 27rpj(A). Below, we will use a special 
notation for the value of the dressed charge on the Fermi boundary Z = Z(±q). A formal 
expression for Z can be given in terms of the resolvent to the operator I — ^=K acting on the 
interval [—q, q] 

9 

Z = 1 + J R(X,±q)dX, (3.4) 



where 



R(\ 0-^Jk(X- /jl)R((i, dfi = ^K(X - 0- (3.5) 

-9 

Finally, we give the formulae for the constants appearing in (|1.5|> . namely the average density 
D, the Fermi momentum k F , and the velocity of the sound on the Fermi boundary vq: 

r 

D = J Pt (X)dX, k F =irD, v o = jf, (3.6) 

-9 

where we denoted e' Q = s' (q). 
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3.2 The dressed energy at low temperature 

We now describe the power-law corrections to the Yang- Yang equation (|2,ip . which appear at 
low but non-zero temperature. For T small enough the solution to the equation (|2.ip e(A) has 
the following expansion |20[ [2T| [22] : 

2 

e(\) = Y,T k e k (\) + 0(T 3 ). (3.7) 

fc=0 

The functions £i(A) and £2 (A) can be found from the analysis of the integral in (|2.ip in small 
vicinities of the Fermi boundaries ±q. The details of this analysis are given in appendix lAl (see 
also [IB [ZD [22]). 

Using (|A.8p and substituting the expansion (|3.7p into the Yang- Yang equation we obtain 



2 2 q 

E Tk£ ^ ~ ^ E Tk f K ( X - aO £ *(m) d l i = \ 2 -h 



k ' - % W* - «> + «(* + «)) - - + op-). (,s) 

What follows from this analysis is that £i(A) = and 

e 2 (X) = -^r{R(X,q) + R(X,-q)), (3.9) 

0£q 



where i?(A, //) corresponds to the resolvent defined by ([3.5 
3.3 The poles of the Fermi weight 



e(A) 

We now consider the low-temperature behavior of the roots r k to the equation 1 + e t = 0, 
namely the solutions to e(A) = 2TriT(k + 1/2), fcsZ. Obviously all these roots collapse to q or 
— q when T — > 0. Therefore setting r^r = g + Tf J + 0(T 2 ) (resp. = —q + Tf^ + 0(T 2 )) and 
substituting these expansions into s(rt) = 2iriT{k + 1/2) we find 



e(q + Tr± + 0(T 2 )) = T^e' + 0(T 2 ) = 2mT{k + 1/2), 
e(- 9 + Tr% + 0(T 2 )) = -Tf±e' + 0(T 2 ) = 2vriT(A: + 1/2). 

Thus, in the linear approximation in T, we obtain two series of roots 

r+ = ± q + ^{k+ 1/2) + 0(T 2 ), A;>0, 



(3.10) 



r 



27TlT ,. , ,,,, 



(3.11) 



A 



±q+^^{k + 1/2) +Q(T Z ), k<0. 



-o 

We will refer to the roots collapsing to +q as the right series and the roots collapsing to —q as 
the left series. 
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3.4 Low-temperature limit of the u(X) integral equation 

From now on, we focus on a fixed contour Cj and consider the associated contribution to the 
asymptotic behavior of the generating function. Therefore, below, we will omit the subscript i 
in the notations of the contour Cj and of the functions Uj(A), Zi(X), etc associated with it. 

Let C be the contour bypassing n points s + in the upper half-plane from above and n points 
s~ in the lower half-plane from below, where n is an arbitrary, but fixed non-negative integer. 
These points s ± are roots of the equation 1 + e~ u<yS ^ T = 0. It is important for our purpose to 
fix the limits of these roots at T = 0. 

The T — > limit of equation (|2.4p coincides with the one of the Yang- Yang equation, hence 
u(X)\ T=0 = £o(A). Therefore it is reasonable to expect that, similarly to the points rt, the roots 

collapse to q or to —q in the T — > limit. Thus, in the low-temperature limit these roots 
should form two series. There will be 

• rip roots s + and n\ roots s~ belonging to the right series; 

• n~ roots s + and roots s~ belonging to the left series. 

Obviously, there exists an integer £, —n < £ < n, such that the numbers rip and are related 
by 

n t + h p = n h + n h =n > n t - n h = n h ~ n v =t ■ ( 3 - 12 ) 
Therefore, for T small enough, one deals with the following structure for the distribution of roots 



s ± : 



{s+} = {q + iTf,+} n+p U {-q + *T^} n -, R(^=) > 0, (3.13) 
{s-} = {q- iTi+} n + U {-q - iT^} > 0. (3.14) 

h h 

The parameters 77, and ^, admit the Taylor expansions fj, = r/. + O(T) and ^, = £, + O(T). 
They appear in the non-linear integral equation (|2.4p defining u(A) and should be computed by 
solving the conditions 

exp (-«*±?&) = ex P (- ^g^) ) = -1. ,3,5) 

Substituting the parameterizations f|3. 13[) . (|3.14p into the equation (|2.4p and expanding up 
to the second order in T we are led to 

u(A) = A 2 - h - £ J K(X-fi) log (l + e" 2 ^ 1 ) d/z + TG±(X) + T 2 G 2 (A) + 0(T 3 ). (3.16) 



Here 



Gi(A) = -2vria - i^(A - /i) d/i, (3.17) 
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and 

(n+ n+ \ / rip nl \ 

j'=i i=i / \i=i i=i / 

It is natural to expect that the solution to (|3.16p has a form similar to (|3.7p 

2 

M (A) = ^T^(A) + 0(T 3 ), (3.19) 

fc=0 

where, as we have already argued, uq(X) = £o(A). Substituting (|3.19p into (|3.16p and using 
()A.10p . we obtain linear integral equations satisfied by the unknown functions lii(A) and U2(X): 

ui(A) - f K(\ — a»)ui(m) ^ = Gi(A), (3.20) 



«2(A) 



1 

2^ 



^(A-/i)u 2 (M)^ = G ! 2 (A) 



7T 



12e; 



r (ET(A-g) + K(A + g)) 



#(A-g)u?(g) A"(A + g )«?(-g) 



4vre' 



4vre' 



(3.21) 



It is then easy to see that the function iti(A) — 27ri^ satisfies (up to a multiplicative factor) 
the equation (|3.3p for the total density at T = 0. As we have already mentioned, on has that 
2-Kpt{\) = Z(X) in the case of the QNLS model, with Z{X) being the dressed charge. Hence, 

Ul (A) = ui(-A) = -2ma t Z(\) + 2m£, a e = a + £. (3.22) 

The solution to equation (|3.21|) can be expressed in terms of the resolvent R(X, fj) (|3.5|1 



u 2 (X) = R(X,< 



2.J>+ + 2vr£e+-^ 

i=i j=i 



7T 



+ uj(q) 



+ R(X,-q) 



n p n h ^ / 2 \ 

.7=1 .7=1 V ^ J 



(3.23) 



It remains to fix the leading Taylor coefficients rj^ and gt. These can be parameterized by 



sets of integers, exactly as it was the case for the roots rt (|3.1ip . More precisely, one has 



u(±q + iTf)±) = ±2mT(p± - i), 
u(±q-iT^) = ^27TiT(hf-^, 



(3.24) 
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(3.25) 



where ui(X) is given by ([3. 221) . 

Remark. Let U\ = ui(±q) = 27ri(£ — a^Z). From now on we assume that u\ satisfies the 
constraint 



Note that the generating function (fl~3j) is periodic over a [23]: ( e 2 ™ Sa; )T = ( e 27Ti( - a+lS > Qx ) T . 
Due to this property the condition above always can be satisfied by appropriate choice of the 
parameter a. Therefore the constraint (|3.26p does not imply any additional restrictions for the 
parameters of the model. 

We stress that the condition (|3.26p is purely technical. It allows us to simplify some interme- 
diate calculations. In particular, it follows from (|3.26p that all the integers and h/j: in f|3.25j) 
are positive. However, one can proceed further without use of the inequality (|3.26p . 

Thus, in this way, we have found the solution u(X) to the equation ([2.4|) up to 0(T 2 ) terms 
and the roots up to OiT) terms. There is no fundamental obstacle to finding higher order 
corrections to u(X) and s^. However, for our purposes, the results obtained here are already 
sufficient. 

4 Correlation lengths 

In this section we compute the rate p[u] of the correlation function exponential decay. In the 
case of the QNLS model, the conformal dimensions giving rise to the critical exponents in the 
asymptotic expansion (|1.5|) were calculated in |24| 125], We now obtain these results by taking 
the T -> limit of equation (f2~T0]) . 

We have already shown in the work [T] how the trivial constant term in (|l,5p can be deduced 
from our approach to the asymptotics at finite temperature. More precisely, this constant stems 
from the contribution of the contour C = R, in other words the case where the sets of the roots 
{s ± } are empty [n = 0). Therefore, in the following, we will only consider the case of non-empty 
sets {s*} (although the results of our analysis remain valid for n = as well). 

By moving the contour C to the real axis, equation (|2,10p boils down to 



The integral over R can be estimated to the leading order in T with the help of (|A.8|) . (IA.10D . 
In its turn, the finite sum is estimated directly by inserting the Taylor expansions of the roots 



7T < 9(ui) < 7T. 



(3.26) 




(4.1) 
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s k . Ultimately, one gets that, to the linear order in T, 



p[u] = -2ia f k F - TM. +tz i J2vt +E^ 7 +£tf 1 + °( r2 )- ( 4 - 2 ) 

o \ j=1 i=1 j=1 j=1 



where we have used (j3.4|) and ()3.6[) . We recall also that u\ = u\{q) = 2ni{l — a^Z). Finally, it 
remains to use that rft, %t are given by (|3,25p . This leads to 



2vrT 

p[u\ = —2iagk F + 



'hp 



(a e Z) 2 -£ 2 -n + Y,Pt + Y,Pj + Y, h t + Y, h J 



3 

j=l j=l j=l j=l 



+ 0(T 2 ). (4.3) 



5 Constant amplitude 

In this section, we compute the low-temperature limit of the constant coefficients Bd[u] (|2.1ip 
and B s [u] (|2,14p . We prove that in this limit, when properly normalized in the temperature, 
B[u] goes to the amplitude of a critical form factor. The latter form factors correspond to 
expectation values of local operators taken between the ground state and excited states where 
all rapidities of the particles and holes are located on the Fermi boundary. We first recall several 
definitions and results concerning the form factors in the QNLS model. The reader can find a 
more detailed exposition in [17jj. 

The form factors of the QNLS model can be parameterized by the rapidities of particles and 
holes [261 |2T [ [2"5 | 129] . If, in the thermodynamic limit (L — > oo) all the rapidities are located on 
the Fermi boundaries ±q, then the corresponding form factor is called critical form factor |17| . 
Hereby the distribution of the rapidities between +q and — q is important. 

Consider a critical form factor parameterized by the rapidities of n particles and n holes. 
Assume that, in the thermodynamic limit, there is (resp. n^) rapidities of the particles (resp. 
holes) going to +q and n~ (resp. n^) rapidities of the particles (resp. holes) going to —q. We 
say that a given form factor belongs to the class, if the numbers h satisfy the conditions 
gathered in (|3.12p . with t being some fixed integer. 

The critical form factors can be presented as a product of a smooth and a discrete part (see 
|17j). The smooth part has a well defined thermodynamic limit L — > oo. The discrete part, 
strictly speaking, has no thermodynamic limit. First of all, it scales to zero as some negative 
power of L, when L — > oo. Second, it not only depends on the rapidities of the particles and 
holes (which are equal to ±q), but also on the quantum numbers associated with the excited 
state. 

In the following, we show that the factor Bd[u] in (|2,lip exactly reproduces the discrete part 
of the critical form factor of the class, provided the distribution ([3. 12|) is fixed. Hereby the 



1 Formally the work |17| deals with form factors of the XXZ spin chain, however the results obtained there 
can be easily reduced to the case of the QNLS model. 



12 



role of large L is played by the inverse temperature: Vo/iT -R- L. The integers and (see 
(|3.24p ) play the role of the quantum numbers describing particles and holes. 

The coefficient B s [u] (|2,14p gives the smooth part of the critical form factor. We first focus 
on the analysis related with B s [u] as the computation of its T — > limit is simpler then for 
B d [u]. 

5.1 Smooth part 

The coefficient B s [u] can be seen as mostly depending on integrals of the following type: 

//= / f'(X)z(X)dX, (5.1) 



c 

where z(X) is given by (|2.9p and /(A) is holomorphic in some domain containing C and M. Then 
moving C to M we obtain 

//-»• j f'(X)z(X)dX-t(f(q)-f(-q)), T^O, (5.2) 

since all roots {s*} go to ±q at T — > 0. Using that, at T = z(A) = for |A| > q and 
z(A) = ui(\)/2iri for |A| < g, we find 

<? 

Urn / /'(A>(A) dX = ± f /'(A)«i(A) dA, (5.3) 



and hence, due to (|3.22p 



9 

y /'(A)z(A) dA = -a, j f'(X)Z(X) dX. (5.4) 



lim 

T->0 

c -i 



Using this prescription we obtain for the limit of the Cauchy transforms 

lim La\z\(w + ijc) = —aiL\_ Q „i [Z](w + i'yc), 7 = 0, ±1. (5-5) 

T—>0 

Similarly 

i 

lim C W = o? / Z ( A ) Z ^) dA dfji. (5.6) 

T^O J (A - fJ,- IC) 2 



-<-! 



Another type of integrals arises in the Fredholm determinant det^ \l — ^K^ u ^ . Recall that 
this operator acts on the contour C as 

I-h Kiu) ] = fW -hi ^ (u) (A,/x)/(/x)d/x, (5.7) 



c 
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where K^(X,fi) is given by (|2.15p . If /(A) is holomorphic in a domain containing C and R, 
then one can easily see that 

J KM(\,riffa)dt* = J K^(X, f i)f(p,)d f i + 0(T), T^O. (5.8) 

C « 

Since it(A) = £o(A) at T = we conclude that, in the T — > limit, the action of the operator 
I — t^K^ coincides with the one of I — ^K^ £ \ The action of this last operator clearly reduces 
to the interval [—q, q] when T = 0. Thus, 

= det[_ M] [J-^]. (5.9) 

Substituting all these results into (I2.14p we immediately reproduce the smooth part of the 
critical form factor obtained in [17] , We give these rather cumbersome expressions in appendix ICl 



lim deti 



Z7T 



lim det 

T->0 



C 



Z7T 



5.2 Discrete part 



The T — > limit of the factor B^u] (|2.1ip is more involved. In order to compute it, we first 
deform the contour C to the real axis. This provides an alternative expression for Bd[u], that 
was originally obtained in |30j 



B d [u] =exp^ J dXdfA ■ (det 



1 



n s j ~ s k 



n , s — 1 / 

x JJ e 2L R W(S7)-2L H W(S+) ^ Ae -2^(A)| j ^ 



-2-iriz(\) 



X=Sl 



Consider the behavior of the Cauchy determinant in (|5.10p at T — > 0. We have 



(5.10) 



det — 



n s j ~ s k 



j>k 

n 

n w-Kf 
j,k=i 



(5.11) 



Now we should substitute here (f3T3]) . <$JM and Hereby at T -> we can set (s^-s^) 2 



T\2 
k) 



Aq , if the roots belong to the different series. Then we obtain 



lim ( T n ~ e2 det - 1 - 



/ / \ -2£ 2 / / \ 2n 

(—l) n + e ( | ( £ 
1 ] \ttJ \2tt. 



n p n h Ti p n h 

n - ^) 2 n (*/ - ^) 2 n (pj - v- k ? n (v - ^) 2 

x - - - - . 



n n(pj*+^-i) 
i=ifc=i 



n n(p7+^- 1 ) i 



j=ifc=i 
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Consider now the low-temperature behavior of the derivatives d\e 2ntz W at A = sf. We 



have 



-2mz{\) 



_uW_ 

1 + e t 
1 + e t 



(5.13) 



Since 1 + e" 



tt(s±) 



we obtain 



Td x e 



-2iriz(\) 



-u'(sf) 



1 + e" 



e(5 j :) ~" (l j ±) 



(5.14) 



Substituting here ([3.19|) . (|3.22p we arrive at 



lim Td\e 
T->o 



-2niz(\) 



e(, (1 - e - 27ria ^)~ 



if belongs to the right series, 
if belongs to the left series. 



(5.15) 



The estimate of the T — >• behavior of the Cauchy transforms Lr[2](s • ) is a more com- 
plicated problem. It is easy to see that one cannot use (|5.5p in this case. Indeed, on the one 
hand the Cauchy transform L[__ q ^[z](uj) on the interval [—5,(7] has a logarithmic singularity at 
u = ±q. On the other hand, we have seen that all Sj go to ±q at T — > 0, therefore Lk[z](s ■ ) 
should diverge as a multiple of logT in the low-temperature limit. In fact, the equation (|5.5p 
only allows one to access to this divergent part; however it does not give an access to the finite 
part of the T — > behavior of the Cauchy transform. 

Similar problem occurs at studying the T — > behavior of the double integral 



A 



(A-^+) 2 



dX dfi. 



(5.16) 



The above double integral can be reduced to a single one 

z(X)z(p) 



A 



(A - M+) ; 



dX dfi 



J (d fl L w [z](n+)) ■ z(fj,)dfj,. 

fc 



(5.17) 



Thus, in order to estimate this integral at T — > one should know the behavior of L^[z](fj,) on 
the whole real axis including the points ±q. This question is studied in details in appendix iBl 
Here we restrict ourselves to present the results of these computations. 

First of all, we list the leading T — > behavior of the exponents of the Cauchy transforms 
£]R[z](Sj )■ There are four cases to distinguish: 

• s£ belongs to the right or to the left series, 

• belongs to the right or to the left series. 
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Using the parametrization ()3. 13|) . ()3. 14j) and (|3.25p we have 
lim LL,[ Z ](± q+ iT^ -1 f i m-Z^ U_A T( Pk ) 



-<j 



.... , ^c^ff) f«4r } = ex P ( / z M- g »+a\ r (^ ± ft 



-9 

Thus, for a given partition (|3,13p . (|3.14f) of roots into the right and left series, we obtain 

1lv,X ( i 

hm n e 2 ^W(^)- 2L «W( s t) = exp |^ | d/i + 2ivma £ Z 



' 1 TO, TO 1 



where we use the standard hypergeometric type notation for ratios of T-functions: 

v q 

-flTia^-llTibk)- 1 . (5.21) 



a\ , ... , dp 
h , ... , b q 



k=l k=l 



In its turn, the exponent of the double integral (|5,16p exhibits the following leading T — > 
behavior 

(/ @)"^ =« a[ * io> (^)- <»■*> 

The functional Ci[.F] appearing above reads 

2 J X — fi J A — q 

-<? -g 

and G(l,x) = G(l + x)G(l — x), where G(x) is the Barnes function. 

Now we substitute ui/2iri = t — otgZ. Combining (|5.12j) . (|5.15|) . (|5.2U|) . and (I5.22[) and using 

1 

C 1 [a l Z{\) -£} = C l [a t Z{\)\ - Ua t [ Z dfi, (5.24) 

-q 

we obtain 

2ra 



Bf ] [u]= lira B d [u] 



g £ 



x Rfd,nt(iP + h{h + }\a t Z - I) R n - n -({p-},{h-}\e ~ a e Z), (5.25) 
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where 



u /( T ruw ^ i >k i >k r 2 \Pk + v), {h k -v) . 

#n,m({p}, {h}\v) = T r t r L t |- ( 5 - 26 ) 

nnfe+^-i) 2 v 

3=1 k=l 



{Pk} , {hk} 



We have thus reproduced the discrete part of form factors given in [17J up to the replacement 
v /T ^ %L. 



6 Final result 

We have calculated the low-temperature limits of the exponential decay and constant coeffi- 
cient in the long-distance asymptotic expansion of the generating function l^e 2ma ^ x )j- for given 
distribution of the roots s^K In this section we sum up all the obtained contributions for I fixed. 

Let us first summarize the results obtained in the previous sections. The asymptotic expan- 
sion of the generating function (e 2l7ra Q-^ T at low temperature has the form 

( e a«raa,) T ~^Q e gto*tk F x x^oo, t -> 0, xT -> oo , (6.1) 
where the coefficients Qi are 

\1 £ oJ 

Here the smooth part £?s [a^Z] is given by (jC.ip . the functional C\ by (|5.23p . To describe the 
factors W± we first define a function W(v, r) depending on complex v and integer r as 

oo ran' 

wm= EE E n^^IK 3 ^ 

™'=o l<pi<— <p n <oo l<hi<— <h n /<ooj=l fc=l 



n,n' : 
n. — n , =r 



2n' 

. SHITTY \ , , , , , 

X(— — ) Rn, n >({p},{h}\v), (6.3) 



where R nn /({p},{h}\v) is given by (|5.26|) . Then 

W± = W(y,r) u=± , z _^ , (6.4) 



r=±t 



It is easy to see that the factors W± correspond to the sums over all the possible choices of 

integers {p^ }™ p and {h^ which parameterize the possible configurations of the roots at 
£ fixed. Such sums are computed in |31] as, 

w< s G\l + r + v) e -^ r(r - 1} 

WM = &oTu) 7 — ^R' (6 - 5) 

v 7 1 — e "o 1 
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Setting here v = ±(agZ — £), r = ±1 and substituting (|6.5p into (J672]) we obtain after simple 
algebra 



2a}Z 2 



^sinh^y (2qZ) 2a t z 
where we have used e' = Zvq. Thus we arrive at the following asymptotic expansion 



{e 2i,a Qx)T „ J- f \ Ai . (6.7) 

l& \ Sm vo J 

Note that the obtained result is obviously a periodic function of a as it was expected thanks 
to the fact that the coefficient Ag is also a function of ag only, see appendix It is also 
straightforward to see that the combination Ag (ttT/vq) 2 ^ 22 coincides with the amplitude of 
the critical form factor of the operator e 2lna ^ x |17| corresponding to the umklapp-type excited 
state of the momentum 2agk F , where Vq/T plays the role of the system size. 

Finally, in order to obtain the long-distance asymptotic expansion of the density-density 
correlation function it is enough to apply the differential operator (|l,4p to the equation (|6.7p . 
Hereby one should distinguish two cases: I = and I ^ 0. In the last case one has due to (jC.ip 



2a? Z 2 



B^[a e Z] =0, d a Bf\agZ] 

a=0 



= 0. 

a=0 



Therefore the second a-derivative should be applied to the coefficient [agZ] . 

On the contrary B^lagZ] = 1 at a = and I = 0, that is at ag = 0, (see [T]). Therefore in 
that case, the second a-derivative should be applied on the combination containing the depen- 
dence on x, namely on e~x.p{2iaik F x] (sinh(7rTx/t>o)) 2oLlZ , as otherwise the second x-derivative 



vanishes. Thus, taking the second a and x derivatives of (|6.7p as it is explained above, and 
neglecting higher order corrections over T as well as subdominant exponential decays in x, we 
arrive at (|l,5p with 

t D 2 £ 2 8 2 ~ 

2 oa z a=o 

It is readily checked that in the x — > oo, xT — > limit equation (|1.5|) does reproduce the long- 
distance asymptotic expansion of the density-density correlation function of the one dimensional 
Bose gas at T = 0, together with the correct values of the amplitudes [19j. 
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A Estimates of integrals with regular functions 

In this appendix we estimate the class of integrals that appears in equations (|2.ip and (|2.4I) . We 
first focus on the integrals of the form: 



J[e] =T J /(A) log (l + e-^) dX. (A.l) 



For our purpose, it is enough to consider the case when /(A) is bounded on the real axis and 
differentiable in vicinities of ±q, although the result of the analysis remains valid at much less 
restrictive assumptions. Due to the properties of e(A) (|3.ip it is clear that 

hmJ[e] = -J /(A)e(A) dX. (A.2) 
-g 

In order to find power-law corrections to the equation (|A.2p . one should estimate the contribu- 
tions coming from the vicinities of ztq more thoroughly. Let 5 > be such that 5 — > as T — > 0, 
while S/T — > oo as T — > 0. We can split the integral J into five parts J = J_ + J- q + Jq + J q + J+. 
The integrals J± correspond to the domains A > q + 5 and A < —q — 5. They behave as 
Q^ e -e{±q±S)/T^ anc j h ence produce 0(T°°) contributions. The integral Jo runs along the do- 

6(A) 

main —q + 5 < X < q — 5. By factoring out e t from the logarithm, we get that 

q-5 

Me] = - J /(A)e(A) dX + 0(T°°). (A.3) 

-q+S 

Finally the integrals J± q correspond to the domains ±q — 5 < X < ±q + 5 and generate all 
power-law corrections in T to (1A.2|) . We now derive the leading power-law correction to (|A.2|) 
coming from the 5- vicinity of q. For doing this, we can replace the functions entering the integral 
by the leading non-vanishing terms of their Taylor expansions. Namely, we replace /(A) by f(q) 
and e(X) by (A — q)e' + Tei(q). Recall that Eo(q) = and we denote e' = £'o(q)- Thence, 

<5 

Jg[e] = Tf(q) J log ^1 + e -^- £l(q ^ dX + h.o.c, (A.4) 
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where h.o.c. means the higher order corrections in T. After changing of variables A = fJ-T/e' we 
obtain 



T 2 f(q) 



Se'JT 



-8e' /T 



log (l + e -^iWj + (fx + £l (q))G(-/i - £l {q))\ d/x 

£1(9) 



(n + ei(q)) dfi + h.o.c, (A. 5) 



where 0(A) is the Heaviside step-function. Using now that 5/T — > 00 we arrive at 



J q [e] 



T 2 f(q) 
2e' 



1 \ fe u 



i / 



00 

J [log (1 + e -») + ^e(-M)] 



T 2 f(q) ( ( , KWfT^m 



Similarly one has 



J-M 



T 2 f{-q) 
2e' 



5e' r 



+ 



6e' n 



+ h.o.c. . 



(A.7) 



Combining (|A.6p . (1A.7D with (|A.3P we obtain after simple algebra 



/WjW//^^, (A . 8) 



2e' 



2e' 



In a similar way, one can obtain the low-temperature expansion of integrals involving the 
function u(X) 

J[u] = T j /(A) log + e"^) dX. (A.9) 

K 

Since uq(X) = £o(A), exactly the same considerations lead us to the estimate 

q 

J[u]=- J f(X)u(X)dX+^(f(q) + f(-q)) 



+ 



T 2 f(q)u 2 (q) + T 2 f{-q)u 2 (-g) + ^ 



2e> 



2e' 



(A.10) 
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B Estimates of integrals with singular functions 
B.l The Cauchy transform in the vicinities of ±q 

In this section we determine the leading T — > behavior of L^[z](A). Its depends on where A is 
located. Recall that 

L RN (A) = ^ AogfUl^ J*!* (B.l) 



2ttz J \ii e -W ) n - A 

R 7 

If A is separated from ±q, then obviously 

where ui(A) is given by f|3.22 j) . 

Let now A — > q as T — > 0. We denote A = X±, if A approaches g from the upper (resp. 
lower) half-plane. Let again 5 > be such that 5 — > as T — > 0, while <5/T — > oo in the 
T — > limit. Consider the contributions to the integral (jB.ip coming from different intervals 
of integration. Obviously when T — > the integrals over domains A > q + 5 and A < —q — 5 
produce exponentially small corrections. On the other hand z(A) can be approximated by 
on the interval [—q — 5,q — 5]: 



2ni 



— 7 log ( 1 + e "?! l ^ — 7' "^l^ , r -> o- (B-3) 

2TTi J + J H ~ A± 2m 7 /x - A± 1 ; 

Extracting the divergent part we obtain 

1 q ~[ 5 u x {n)dn ^ 1 fttiW-1.1 ^+^ T logf A± - g + ^, 5^0, (B.4) 



27TZ y /X — A± 27TZ y /i — A± 27TZ \ A± + 

and we remind that U\ = u\(±q). 

It remains to compute the integral over [q — 5,q + 5]. Following the method of the previous 
section we linearize the functions u(/j.) and s(/j.) in the vicinity of /i = q. Then we have 

Replacing /x£q/T = £ we arrive at 



/g= 2^ 

-fc^/T 



+ ^Z7 lo S jcr h ( B - 6 ) 



e-t± 2vn b w ± + ^ 



T 
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where c is an arbitrary positive constant and we have set t± = (A± — q)e' /T. We can send now 
5/T —7- oo. Substituting into (1B.6|) 



±i / e =F<"«-*±) dw, 



(B.7) 



we arrive at 



oo oo 



2tt 



dl; J dio 

oo 



i og (^-+jL^)+ Ul e(-£-c) 



3 =FM£-*±) 



+ 2^ l0g 



t± + c 



(B.8) 



The integral over £ can be calculated by means of an integration by parts followed by a compu- 
tation of the residues at + e~ Ul = and + 1 = 0: 



log 



e« + 1 



+ Ui0(-£-c) 



1 - e ±iuUl u v r 



2u sinh(7Rj) 2-kiuj 



(B.9) 



Thus, we arrive at 



T Ul 1 

^ = 2^ l0g 



t± + c 



2ir I uj 



7T 



sinh(7rw^ 



(1 



(B.10) 



Due to (|3.26p the last integral is convergent. It can be computed in terms of the T-functions via 



e'P" duj 



b — a 



TT 



sinh(7ro;) 



e' aui - e- bu) 



(a—b) log 



2tt 



fp+b | 1 
^ 2tt "I » 



+^og } . (B.ll 

1 \ 2tt ~1~ 2) 



Thus, we obtain 



t Ul l 
^ = 2^ l0g 



±2iriT 



(X±-q + S)e' 



Y i l _|_ (^±~g) e o _|_ jui. 



7 ± log 



2mT 



2-i 



± 



(A±-(?)£p 
2?riT 



Combining this result with (|B.4p we find the following estimate 



(B.12) 



Lr[z](X±) = — 



log 

2vri 



- «l , «i, ( X±-q 
dji + — log 



(i-\± 

(A± - g)gp 
±2vriT 



27ri \ A± + g 

1 i (A±-g)£p 



±log- 



1 1 2 ^ 27rir ^ 27TJ 



pf 1 i (A±-<?)£q 
M 2 ^ 2?riT 



T^O, A~g. (B.13) 
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Similarly, if A— > — q as T — > one has 

g 

j [ i/> \ 1 fui{n)-ui m {X±-q\ 

Lr[z]{\±) = — / d(i + — log - — — 

2m J [J> — a± 2m \X±+qJ 

-q 

, , . -n ( 1 i (A±+g)£p \ 

^1 . f(\± + q)e? Q \ l \2± 2-KiT +2ni) 



B.2 The double integral 

Consider now the low temperature behavior of the integral A given in (|5.17p . As usual we split 
the integration domain into several pieces A = A^ q + Aq + A q + A + . The integral A + (resp. 
A-) over the domain A > q + 5 (resp. A < —q — 5) are again exponentially small in respect to 
the T — > limit. When A 6 [—q — 5,q — 6], we can use the expression (|B.2p for Lr[z](A+) and 
also replace z(X) by ^ ) . This gives 

q-S q 

A « s ^ I < B - 15 > 

-<7+5 -g 

Integrating by parts we arrive at 

Ao = _i 7 5 dA ? d m^M + ^zJl ( f M»)d» Jm>^\ 

(2m) 1 J J A*-A + (2my \J - g + <5 7 /x_ + g - i5 I 

-g+<5 -g Vg -i? / 

Here we have used that ui(X) = ui(— A). This last property also allows one to symmetrize the 
integrand, so that upon sending 5 — > 0, we get 

7* dx I d ^m^t> = ' / "iM"iM - ..OHM ^ 

(27rz) 2 7 7 M - ^+ 2(27rz) 2 7 A - /i 

-g+i5 -g -g 

(B.17) 

Extracting then the divergent part from the single integrals in (|B.16P we find 



MX) 

2m 



2u 2 ( S \ 



where the functional Ci[-F] is defined in f|5.23j) . 

Consider now the contribution A q coming from the interval q — 5<X<q + 5: 

q+S 

A q = J (d x L R [z](X + )) -z(X)dX. (B.19) 

q-S 
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Substituting here ()B,13P we arrive at A q = + A K q l) + Af, where 

q+5 

z(X)d x 



,(2) 



,(3) 



A q 



1 

2tH 



q-5 



mdi) -u x u\ . . 

x dp--?—, log A± + q) 

fj, — A± 27H 



dA, 



(B.20) 



and 



q+8 



q " (2vri) 2 



log 



-q-5 



"(*) 

1 + e t 
1 + e t 



dA 



(B.21) 



4 3) = — 
9 2vri 



q+S 



d\ < lop 



1 , (A+-(?)£q 

2 "i" 27riT 



+ 



2-i 



-q-5 



p / 1 , (A+-i?)£b 
1 1 2 " r 2vriT 



"1 

2vri 



log 



(A+ - q)e' 



2mT 



x log 



"(A) 

1 + e t 
1 + e r 



dA. (B.22) 



It is easy to see that A^ -)• as 5 -)• 0, because the integrand is a bounded function as 
T — t- 0. The integral .Ag can be estimated similarly to (|B.5P : 



u 



A (2) _ 

q (2vri) 2 



^log 



-2mT 



log' 



'I _ JUJ 
,2 2tti ) 



5e' J 2ni 6 r(|) ' 



0, 8-+0. 



(B.23) 



As for the remaining integral A q 3 \ its leading behavior is obtained by a linearization of the 
functions u(X) and e(A) in the vicinity of A = q. After the change of variables £ = (A — q)e' /T 
followed by an integration by parts, we find in the 5 — > 0, T — > limit 



4 3) = — 

q 2m 




nil £±«i 

1 1 2 "1" 27TJ 



2 T 27TJ 



«1 , 



2vri 




1 



1 



1 + e-t~ Ul 1 + 



d£. (B.24) 



We close the integration contour in the upper half-plane and compute the integral (|B,24p by 
residues. These are located at £ = —U\ + iri(2k + 1) and £ = iri(2k + 1), k = 0, 1 . . . . Hence, 



4 3) - E 



fc=i 



log 



r 2 (fe) 



r(fc + ^)r(fc-^) 27ri 

V 2m/ \ 2-Ki! 



Ul 

7T-l°g 



b — I — "1 

_ 2 27Tt 

fc- 5 



+ lim 



log 



G 2 (N + 1) 



^l log r (^-^) 



(B.25) 
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where G{x) is the Barnes function and G(l,x) = G(l + x)G(l — x). Using the asymptotic 
behavior of the T and Barnes functions for z — > oo with z R_ 



logG{z + 1 + a) - logG(z + 1) = alogV27r + \{2z + a) log z - az + o(l) 
log r(z + 1 + a) - log T(z + 1) = a log z + o(l), 



we find that the limit in the last line of (|B,25P vanishes. Hence, 



Combining this result with (|B.23p we find 



log 



-2mT\ , / m \ 



T->0, 5->0. 



(B.26) 



(B.27) 



Similar calculation in the vicinity of the point — q leads us to the following below contribution 



coming from the interval q — 5 < X < q + 5: 

-q+8 

A_ n = 



j (^H(A + ))-,(A)dA^^lo g (^) + lo gG (l,|l), T^O 

-q-S 



Thus, taking into account (|B.27p . (|B.28|) and (|B.18|) we finally obtain 



Ci 



m(A) 

2vrz 



'(is)* gj— 



_U]_\ 



0. 

(B.28) 
(B.29) 



C Smooth part of the amplitude 



In this section we give the exact expression for the smooth part of the amplitude Bg[aiZ] 
]hn.T->o B s [u]. Provided the condition (|3.12|) holds, we have 



Bf\ ai Z] = {e 



2nia 



2 o -C . 



det ( 1 + 2¥?^ (1) det ( 1 + 2^ (2) 



( det [I-±K]Y 



x ^ e -a e L[Z](9i+ic) _ e 27ria-a £ L[Z](6»i-ic)j ^ e a t L[Z](0 2 -ic) _ e 2nia+aiL[Z\(e 2 +ic)\ _ ^Qjj 

Here L[Z](w) stands the Cauchy transform of the dressed charge Z(\) on the interval [— q, q], 
and Co is given by (|5.6p . The integral operator I— -£=K acts on the interval [—q, q] and its kernel 
was defined by (|2,3p . The operators I + -^U^ [a^Z] and I+-^U^[agZ] act on a anticlockwise 
oriented closed contour surrounding [—(7,(7]. Their kernels are 

f^W, M) = -e-^^H . K ^- W ' ) - 2 Kai9l - z W ' ) y (C.2) 

v ' p—aiL[Z\(w+ic) g27r«a— aiL\Z\{w— ic) 1 y ' 
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and 

#< a w, m) = • f t (ro ; m,) ^; (ro J 2 t, (c.3) 

v ' » L «• j/ e aiL[Z](w'— ic) _ e 2Tna+aiL[Z](w'+ic) ' v ' 

where X a (A) is given by (|2.18jl . Finally parameters 9\ and #2 are arbitrary complex numbers 
lying inside of the contour where the operators CA 1 ' 2 ) (w,w' , [cegZ]) act. If we set 6\ = —q and 
02 = q, then we reproduce the smooth part of form factors of the class obtained |17j . 
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